Nonequilibrium Spin Dynamics in the Ferromagnetic Kondo Model 
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Motivated by recent experiments on molecular quantum dots we investigate the relaxation of 
pure spin states when coupled to metallic leads. Under suitable conditions these systems are well 
described by a ferromagnetic Kondo model. Using two recently developed theoretical approaches, 
the time-dependent numerical renormalization group and an extended flow equation method, we 
calculate the real-time evolution of a Kondo spin into its partially screened steady state. We 
obtain exact analytical results which agree well with numerical implementations of both methods. 
Analytical expressions for the steady state magnetization and the dependence of the long-time 
relaxation on microscopic parameters are established. We find the long-time relaxation process to 
be much faster in the regime of anisotropic Kondo couplings. The steady state magnetization is 
found to deviate significantly from its thermal equilibrium value. 

PACS numbers: 72.10.Fk 72.15.Qm 72.25.Rb 73.63.Kv 75.20.Hr 75.30.Hx 76.20.-|-q 



Introduction. Recently it has become experimentally 
feasible to trap isolated single molecules in nanogaps 
forming transistor geometries. In such molecular quan- 
tum dots a variety of interesting new phenomena have 
been observed. In case of a single Cgo-molecule attached 
to metallic leads (sketched in Fig. 1) the quantum phase 
transition between a singlet and a triplet eigenstate of 
the molecule has been studied in detail [1, 2, 3]. In par- 
ticular, if the isolated molecule is prepared in the triplet 
configuration, its spin is partially screened by the conduc- 
tion band. In this case the resulting effective exchange 
interaction between the residual spin and the conduction 
band is known to be ferromagnetic [4, 5]. 

Replacing Ceo by a single- molecule magnet (SMM) 
such as Mni2 gives rise to even more complex quan- 
tum impurity physics [6]. As a result of magnetic 
anisotropy induced by spin-orbit coupling, the large in- 
trinsic spin of the SMM tends to align along the easy 
axis of the molecule. This gives rise to an energy bar- 
rier which suppresses magnetization reversal and makes 
SMMs promising candidates for applications such as 




FIG. f : Sketch of a Ceo-molecule coupled to metallic leads. 



high-density magnetic storage and quantum information 
processing [7]. When coupled to metallic leads, the SMM 
can be described by an effective Kondo Hamiltonian with 
anisotropic exchange coupling between the impurity spin 
and the conduction band [8, 9] 
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In the cotunneling regime it has been shown [10] that the 
exchange interaction is ferromagnetic, i. e. Jy <§C Jl < 0, 
if adding or subtracting an electron to the molecule in- 
creases the spin of the SMM. Preparing the system in a 
well-defined spin state and measuring the real-time spin 
dynamics can be achieved using electrical [11] or optical 
[12, 13] field pulses, albeit experimental challenges in ap- 
plying these techniques to molecular quantum dots still 
remain. 

Although the Kondo Hamiltonian has been studied in 
great detail, the ferromagnetic regime has often been ne- 
glected (an exception is the investigation of the spatial 
equal-time spin-correlations of an underscreened spin-1 
impurity [14]). In this letter we will focus on two impor- 
tant questions arising in this context: By studying the 
magnetization dynamics we will investigate how fast an 
initially polarized spin will reduce its magnetization due 
to spin flip scattering. For the antiferromagnetic Kondo 
model this question has been answered in [15, 16]. Our 
analysis yields important information about the domi- 
nant relaxation mechanism in related experiments. 

One further important question regarding the relax- 
ation process is the nature of the final state of the quan- 
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turn system. Since a pure state remains pure under uni- 
tary time evolution, the complete system is not expected 
to behave like an equilibrium state even at long times. In 
equilibrium, the conduction electrons will not fully screen 
the spin, as is well known for the ferromagnetic Kondo 
model [17]. Instead, the coupling Jj^ provides weak spin 
flip scattering which renormalizes the magnetization of 
the impurity spin to some finite value in presence of a 
symmetry-breaking infinitesimal magnetic field. In the 
case of isotropic couplings and a spin S, this value is 
known as [17] 

{S.)=S{l + ^ + 0{j')), (2) 

where p is the density of states in the conduction band 
with support [—D,D] (we assumed p — {2D)^^ and em- 
ploy units in which h = kg = D = 1 in the following). 
However, the fact that the system's final state differs from 
thermal equilibrium does not necessarily imply that lo- 
cal observables retain a memory of the initial prepara- 
tion, and in fact even for certain integrable systems the 
reduced density matrix of a local subsystem is known to 
thermalize [18]. For the ferromagnetic Kondo model, the 
low-energy spin-fiip scattering rate renormalizes to zero, 
leaving open the question whether in the steady state 
at long times the impurity has finite magnetization or 
not. We will give a definite answer to this question in 
this letter and show that the asymptotic nonequilibrium 
magnetization at long times differs from the equilibrium 
value. Therefore in our model information about the ini- 
tial preparation of the system is never completely lost, 
even for local observables. 

Methods. In recent years, several numerical approaches 
have been developed to calculate real-time dynamics of 
quantum impurity systems [16, 19, 20, 21]. However, the 
accuracy of numerical data is usually not sufficient to 
give precise answers about the nature of the long-time 
decay, i. e. to identify analytical laws for the long-time 
tails and steady state values. We therefore use an ana- 
lytical approach to identify the long-time behavior and 
compare it against numerical calculations to validate our 
analytical approximations. We first describe our ana- 
lytical approach, before briefly sketching the numerical 
technique. 

Within a poor man's scaling analysis a ferromagnetic 
exchange coupling of an impurity spin to a fermionic bath 
renormalizes to zero at the Fermi energy. This allows per- 
turbative renormalization techniques to accurately de- 
scribe the low-energy physics of such a system. In this 
context, a powerful technique is the flow equation method 
as invented by Wegner [22] and independently by Glazek 
and Wilson [23] . In a recent modification of the original 
flow equation method it has been shown that the underly- 
ing renormalization scheme can be extended to calculate 
the real-time evolution of interacting many-body systems 
[24, 25, 26]. As a notable feature, this approach allows to 



derive exact analytical results. For further details of the 
flow equation approach and its application to the Kondo 
problem we refer to Ref. [27]. 

We briefly outline the main steps of the flow equation 
calculation. Details of this calculation will be published 
elsewhere [28]. As usual, the impurity spin operator is 
first transformed by a sequence of infinitesimal unitary 
transformations. The flowing spin operator has the form 

S.{B) = h{B)S.+Y,lk'k{B) : {S+s-,^+S-s+,^) :, (3) 
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where the initial form of the operator is obtained for 
5 = [27]. Here, the operators s^^. are matrix elements 
of the conduction electron spin density raising and lower- 
ing operators. The coupling constants h{B) and "fk'k{B) 
obey the flow equations 

^ = Y.{s,,-e,)Ji^,{Bh,k'{B)n{k'){l-n{k)) 
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^ = hiB){e,, - e,)Jt,,{B) + 0{J'), (4) 

where n{k) denotes the Fermi distribution function. In 

addition, the flowing couplings J^'ki^) ^'^'^ "^I'ki^) 
the Hamiltonian enter, which have to be calculated sep- 
arately [27]. The fixed point of the transformation is 
reached in the limit B oo, where we denote coupling 
constants in this basis by a tilde, e. g. h. The Heisenberg 
equation of motion of the impurity spin can be solved 
efficiently by solving it first for the transformed impurity 
spin and reverting the unitary flow afterwards [24, 25]. 
By solving the equations of motion in the diagonal basis 
of the Hamiltonian, one avoids secular terms that grow in 
an uncontrolled way with time and can obtain controlled 
analytical results even for the asymptotic long-time be- 
havior. 

In order to verify this semi-analytical approach we em- 
ploy the recently introduced time-dependent numerical 
renormalization group method (TD-NRG) by Anders and 
Schiller [16]. Describing this method in detail is beyond 
the scope of this letter. Let us only mention here that 
it is tailored to calculate the response of an (arbitrary) 
quantum impurity system to a sudden quench at time 
t = and is able to access the long time-scales charac- 
teristic for Kondo physics. We refer the interested reader 
to Ref. [16] for more details on this method. 

Results. For our analytical calculations, we assume 
that the impurity spin is prepared in the up state of the 
spin projection operator Sz before the thermalized con- 
duction electron bath ([-FS')) is coupled to it, leading to 
a product initial state 

|^) = |T)®|i^5). (5) 

At time t = 0, the spin is coupled to the conduc- 
tion electrons, which would e. g. be realized by attaching 
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metallic leads to the single molecule magnet. In the fol- 
lowing, we restrict our calculations to spin 5* = 1/2 and 
zero magnetic field. Numerical calculations which we per- 
formed for the dynamics of larger spins S all satisfied the 
trivial relation (S'z(t)) — S{az{t)), which is exact up to 
0{J'^S^) from higher order flow equations [28]. Here S 
denotes the size of the spin and az is the spin operator 
for spin S — 1/2. After solving the Heisenberg equation 
of motion for the operator Sz, the formal result for the 
magnetization reads 
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(6) 

In the following, wc will only discuss the purely quan- 
tum mechanical case, corresponding to T = 0. It turns 
out that the magnetization dynamics can be fully under- 
stood by the energy dependency of the couplings "fkk', 
which are obtained from a solution to Eq. (4). 

Isotropic Kondo Model. Let us first investigate an 
isotropically coupled spin. In equilibrium, perturbative 
scaling shows [29] that the isotropic coupling J logarith- 
mically decreases upon reducing the half band width from 
D to some A < Z?, 
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At the low energy fixed point, an infinitesimal magnetic 
field is sufficient to polarize the free spin, leading to a 
finite magnetization according to Eq. (2) . 

In nonequilibrium, our results show that the magneti- 
zation saturates as well, but to a different value than in 
equilibrium. Using the exact low energy behavior of the 
couplings 7fc/c' in Eq. (6), the asymptotic behavior of the 
magnetization is obtained as 
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This behavior can be understood from the logarithmic 
renormalization of the coupling J, which directly enters 
the low energy flow of the couplings ^kk' via Eq. (4). 

The steady state magnetization {Sz{t oo)) = i(l -|- 
pj ~\- 0{J^)) therefore differs from the equilibrium value 
as given by Eq. (2). The reduction from full polarization 
is pJ/2, which is twice the equilibrium value. This can be 
attributed to the fact that the nonequilibrium dynamics 
starts with an impurity spin that is not dressed with a 
conduction band electron cloud: it therefore relaxes to 
a smaller value of the magnetization as compared to the 
dressed impurity spin in equilibrium. 

A direct numerical solution of the flow equations al- 
lows to accurately determine the relaxation process also 
at intermediate and short time scales. Together with the 
analytical result from Eq. (8), this calculation can be 
compared to TD-NRG calculations. Both methods yield 
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FIG. 2: Results for the isotropic ferromagnetic Kondo model 
(x: TD-NRG data, □: flow equation data). Using our ana- 
lytical result we fitted our data against {Sz{t)) = (1 + aj\\ + 
(ln(t) — l/(cj|| ))^"'^)/2 using a, c as fit parameters (lines). 



very good agreement up to time scales of order i w 10'^ 
where the asymptotic logarithmic relaxation is clearly 
visible, see Fig. 2. Increasing deviation of the curves 
for larger coupling strength J can be explained by the 
O(J^) corrections to the flow equation result, which we 
neglected. We checked that the relative deviation of the 
two methods in terms of the quantity {Sz{t)) — 1/2 at 
some large but flxed time indeed grows approximately 
linearly in J. A flt of the TD-NRG curves and the nu- 
merical implementation of the flow equation approach are 
in good agreement with the analytical result of Eq. (8) . 

Anisotropic Kondo Model. Studying the anisotropic 
Kondo model we restrict ourselves to the experimentally 
relevant case Jj_ > Jy in the following. From a poor 
man's scaling analysis [29], it is known that the coupling 

p /j2_ j2 

Jl renormalizes as Jj_(A) cx A V n ^ at low energies A. 
As in the isotropic case, the flow equation analysis shows 
that this behavior determines the asymptotic long-time 
relaxation of the spin, given explicitly by the power law 



{Sz{t))^ 0.51 1 



2.911 



2511 



0(J2) 



(9) 



where ^|| = —p^J^^—Jj_. The constant a derives from 

the scaling equations for J± and J|| . Numerical checks 
show that it can be replaced by a w pj± as long as 
J|| < 2Jj^. In comparison to the isotropic case, the power 
law decay of spin flip scattering at low energies leads 
to much faster relaxation of the magnetization, whereas 
the steady state magnetization is enhanced. This behav- 
ior is reproduced by our numerical calculations shown in 
Fig. 3. Again, our calculations showed that the steady 
state magnetization {Sz{t — > oo)) = 1/2-1- (a^/(4^||) is 
reduced twice as much from full polarization than in equi- 
librium. The analytical results are confirmed by numer- 
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FIG. 3: For the anisotropic ferromagnetic Kondo model 
our numerical findings coincide with our analytical results 
using both methods. Fitting our data against {Sz{t)) = 

_ /j2_ j2 

a • t V II ^ + c we found good agreement for the fit pa- 
rameters a, c. 



ical fits of our data, see Fig. 3. Let us point out that for 
the anisotropic Kondo model, our methods are starting 
from slightly different initial states. Using the flow equa- 
tion approach one is restricted to a situation where the 
spin is initially completely decoupled from the fermionic 
bath. On the other hand, stability of the TD-NRG al- 
gorithm in the anisotropic model requires preparing the 
polarized spin at time i < by applying a large magnetic 
field, while still allowing for a small exchange coupling to 
the metallic leads. The same long-time power-law relax- 
ation was obtained with both methods. However, the 
slightly different initial states used in both methods be- 
come significant on short and intermediate time scales. 

Conclusions. We employed two different methods to 
analyze the real-time evolution of a ferromagnetically 
coupled Kondo spin, which is initially prepared in a po- 
larized state. Exact analytical results for the long-time 
behavior of the magnetization were obtained for two dif- 
ferent situations. For the isotropic ferromagnetic Kondo 
model, the long-time relaxation is logarithmic in time, 
whereas anisotropic couplings lead to a power-law de- 
cay at large times. Furthermore, exact analytical results 
for the asymptotic nonequilibrium magnetization were 
presented, which differ from the equilibrium magnetiza- 
tions. They confirm that the local quantum impurity 
retains a memory of the initial preparation for asymp- 
totically large times. This is due to the combined effect 
of nonequilibrium preparation and ergodicity breaking 
already in the equilibrium system. 
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